Exactly solvable model of a quantum spin glass 
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A mean field spherical model with random couplings between pairs, quartets, and possibly higher 
multiplets of spins is considered. It has the same critical behavior as the Sherrington-Kirkpatrick 
model. It thus exhibits replica symmetry breaking. The order parameter function is solved exactly 
in the whole low temperature phase. The zero field cooled susceptibility remains finite at low T. 
Next a quantum version of the system is considered. Whereas the magnetic properties are not 
altered qualitatively, the thermodynamics is now regular at small temperatures. 
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The phenomenon of spin glasses has posed a major new field of research. The basic phenomenon is breaking of 
ergodicity. This leads to highly non-trivial properties as observed recently e.g. in aging experiments [1] The most 
known spin glass model was formulated by Sherrington and Kirkpatrick (SK) [2]. Its solution problem was presented 
by Parisi, see [3] for a review. It gave a description of ergodicity breaking in terms of overlaps. For more recent 
reviews on spin glasses, see [4], [5], [6]. 

The spherical model has always played a special role for understanding the basic phenomena of phase transitions. 
It was introduced by the late Marc Kac and solved by Berlin and Kac. [7] It constituted a simple model for studying 
ON critical behavior, see Joyce for a review. [8]. 

The spherical random bond spin glass was studied in ref [9]. It shows freezing but no replica symmetry breaking. 
This model was extended by the present author to include short range ferromagnetic interactions. [10] The spherical 
model with random couplings between p spins was studied recently by Crisanti and Sommers (CS). For p > 2 it was 
found that there occurs a one-step replica symmetry breaking. [11] Dynamical aspects of this transition were studied 
in detail. [12] 

For Ising spin glasses the Parisi solution is known explicitly near the spin glass transition. In the frozen phase 
the solution has been put in the form of a stochastic differential equation, for which no explicit solution is known. 
As a model solvable in the whole low temperature phase we propose here a mean field spherical model with random 
couplings between pairs, quartets, and, possibly, higher multiplets of spins. In contrast to previous spherical models, 
where only one type of these couplings occurs, the model belongs to the universality class of the Ising spin glass. We 
shall derive the explicit form of the order parameter in the low temperature phase. Then we shall consider the recent 
quantum formulation of the spherical model [13] and study its low temperature behavior. 

Consider a system with N spins with random couplings between sets of p spins. The Hamiltonian reads 
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The J's are independent Gaussian random variables with average zero and variances (J^^...^) = (p — ly.JpN 1 p . 
The spins are subject to the spherical constraint ^ . Sf = Na. The classical partition sum reads 



Z d = / D5e-"« (2) 



where DS = Sf — Na) Ylii^^i/V^)- The thermodynamics of this model with only one of the p terms was 

worked out by CS. In a replica formulation of the average partition sum it was found that the free energy only depends 
on the overlap q a p = (l/N) . This result has to be summed over p. We introduce 
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and we obtain a replicated free energy 
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If J2 = the model has the one-step replica symmetry breaking solution studied by CS. When J2 > J* z (J3, J4, • • •) 
is large enough, a continuous transition will take place. From now on we assume that this is the case. In eq. (4) 
also the diagonal terms q aa = q d occur; in the classical situation they are equal to q d = c Expanding in off-diagonal 
elements q a p one obtains 

2 pF n = -n{\ + log q d + p 2 f(q d ) + p 2 H 2 q d } - /3 2 H 2 ]T ' q a p - \(p 2 jf - \) ]T ' ~ 2 
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where the primes indicate exclusion of diagonal terms. In the case J3 = 0, but J4 > 0, this is exactly the relevant part 
of the free energy functional of the SK model [14] [15]. The model thus belongs to the same universality class, and 
exhibits replica symmetry breaking. This continuous transition sets in at Tq = J'lld = J'i&- F° r J3 > the model 
belongs to the universality class of the random bond Potts model. When the J2 > J| and the next non-zero term 
is J p for some p > 4, the model belongs to a new universality class. We shall not consider such situations here and 
assume J 2 > , J3 = 0, J4 > 0. The higher couplings are only relevant well below Tq. 

We express the q a p in the Parisi function q(x). It takes the plateau value q\ for x\ < x < 1. In an external field it 
also has a plateau value q(x) = qo for < x < xq. The inverse function is x(q). Using the expressions of an appendix 
of ref . [11] we obtain the explicit expression for the classical free energy 

2pF cl = -p 2 [ dx{f(q d ) - f(q(x)) + H 2 q d - H 2 q(x)} - - \og(q d - q,) - 1 (6) 
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where 



The saddle point equation for q(x) reads 



I(q) = Qd-qi+ / x(q')dq' (7) 
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P 2 f(q(x))+P 2 H 2 = f (X) dq'-f- (8) 
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In the region where q'(x) ^ one has j3 f" (q) = I(q) . It follows that x(q) has a universal shape at all T 

which gives q(x;T) = q(/3x) after inversion. In an external field H the plateau value qo follows from 

H 2 = q f"(q )-f'(q ) (10) 
This condition is independent of T. q\ follows from 

T 
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Of special interest is the case where the couplings J p are such that 

/(*) = '!flog^ (12) 
4a 1 — aq 

for some a in the range < a < I /a. In zero field one finds q(x) = (j3J'2/'2a 2 )x for < x < x\ and 

qi = (13) 

J2 + y/J% - 4a 2 T(q d J 2 ~ T) 

Note that here the order parameter function q(x) is either linear or constant in the whole low temperature phase. The 
breakpoint x\ = 2a 2 Tqi vanishes both near the critical temperature and near zero temperature. 
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In an external field the freezing transition is given by the AT-line [16]. Here the replicon mode [15] becomes 
massless. We propose to give this mode the more physical name ergodon. The critical field where this occurs, Hf r (T), 
is obtained by inserting qo = Qi(T) into eq. (10). Note that Hj r remains finite at T = 0. 

In the high temperature phase the susceptibility follows the Curie law x = fild- The zero field cooled susceptibility 
reads 

xzFc=p{qd - qi)= 7fkm (14) 

It remains finite at T = 0. This is usually observed experimentally. However, the SK model predicts a vanishing value 
at T = 0. [17] The field cooled susceptibility 

f 1 1 

XFC = P I dx(q d - q(x j) = — (15) 
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is constant in the frozen phase, in good agreement with experiments. The internal energy and the entropy read 

U cl = /3[/( 9 i) - f(q d )] + " f * dq^Fiq) (16) 

2 S cl = f [f( qi ) - f(q d )] + 4== + 1 + l°gfe " <Zi) 
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For small T one has S r j ps^ — \ log/3 2 /"((/i). As usual for classical vector or spherical spins, it goes to — oo as T — ► 0. 
This also implies that the specific heat goes to a constant, C —> 1/2, at low T. For the case a = 1/2 with pair and 
quartet couplings (J 2 = 1, J4 = 2), the specific heat and the entropy have been plotted in Figure 1. 

We have wondered whether the anomalous low temperature behavior of our model can be cured. In order to do 
so, one should regularize the low temperature behavior of the spherical model in general. This can be done by going 
to a description of spherical spins in a standard thermal field theory for bosonic quantum spins. [13] The situation 
becomes simplest when the spins are complex valued and couplings are Hermitean. For the case of random pair and 
quartet couplings we now assume the Hamiltonian 

= "IE J ij% S j - E JijktfSjSlSt - H ]T(ReS 8 + ImS,-) (17) 
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For each pair there are now two independent random variables, J' and J", with average zero and variance 

j|/2#, in terms of which J 8 j = J*- = J' + For each quartet k, I) with i < k, j < I there are four independent 
random variables, J[ 2 and J" 2 , each having average zero and variance 9J|/2# 3 . In terms of J 1<2 = J[ 2 + iJ" 2 the 
couplings in eq. (17) read Jijki = Jjuk = Ji + 1J2 an d Jijik = Jjiki = Ji ~ *'^2- The thermal partition sum involves 
spins on the imaginary-time interval < r < [3 

Z = JvSe^jyr^^SUr^S^-niSAr)^ (18) 

with boundary conditions Si(/3) = Si(0). The r-axis is discretized in M = /3/dr steps, with first N —> 00 and then 
M —> 00, while cIS(t) = S(t + dr) — S(t). The integration measure is a repetition of the spherical measure at all r, 
VS = C N Y\ T DS(t), where DS = 6(2Na — S* Si)Hi(dS* dS{ /2ir) is the spherical measure for complex spins and 
C= (l/2a)n n7£0 (7r|n|/2a). 

For a ferromagnet eq. (18) leads to similar critical behavior as the classical spherical model. However, the entropy 
is non-negative and vanishes at T = 0. The zero point magnetization, Mo = \/ a — a, shows a quantum reduction 
from the classical value Mo = \[g . [13] 

Here we wish to see whether the low temperature behavior is also cured for our spin glass model. We therefore 
extend the Crisanti-Sommers approach to thermal fields. We introduce a "Fermi level" \i a (j) related to the spherical 
constraint at imaginary time r and overlaps q a p{T, r') =(l/N) S* (t)5 , 8 '(t'). We look for time-invariant solutions 
where /j„(t) = /j, and where q a p{j, r') = ^ e q a f}{t) exp{27rinT(r — r')} with bosonic Matsubara frequencies e = 
TrnT/2a. In the high temperature phase at zero field only the q aa (e) are non-zero. We shall denote q d = q aa {e = 0), 
which now becomes a smooth, increasing, strictly positive function of temperature. The spin glass temperature Tq is 
the solution of the relation Tq = J'z^diTa)- Here the off-diagonal elements q a f)(e = 0) become non-zero; since Tq < <r, 
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this occurs below the classical transition point. The q a f}{t) with e ^ are always equal to zero. [18] For f / Owe set 
p c = fiq aa {e) and find that the full free energy reads 

F = 2F d (q(x); q d ) + F qc (q d ; n;p e ) (19) 



where F r j is the free energy given by eq. (6) with f(q) = J 2 q /2 + J 2 q /4. (The factor 2 in eq. (19) 
spins are complex.) The quantum correction equals 



arises since now 



(J 2 ] T 2 J 2 
PF qe = - \ -fpl + !"(/« + 2ie)p e + log 2iep £ [■ ^ ^ P^P^P^Pe, + M«d " + log 2 « (20) 

where terms with all e 8 ' equal to zero are excluded; they occur already in F r j. In the quantum situation qd, p, and 
all p c have to be treated as additional variational parameters. For large temperatures one finds p = Ty/a, qd = a/y, 
p c = l/(/U + 2ie), with tanht/ = a/<7. The entropy Soo = yc / a — log2sinht/ agrees with the result in the situation 
without interactions. [13] Only in the classical limit a — ► one finds back previous relation qd = a. Therefore the 
susceptibility x = Qd(T)/T has been reduced from its classical value at any T. 

First consider the situation where only pair couplings occur, J4 = 0. This case, where no replica symmetry breaking 
occurs, is exactly solvable since only Gaussian integrals occur. The free energy may be expressed as 

f3F = -Pup + / rfA/9(A) log 2 sinh a/3( / u - A) (21) 
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Here p(\) denotes the semi-circular density of eigenvalues of the coupling matrix J 8 j, viz. p(\) = \/4J| — \ 2 /2ttJ 2 
Alternatively we can solve the p e from eq. (20). They read p e = l/[^/U + ie + \J {\p + «e) 2 — J|] with [i = 2J2 in 

the condensed phase. When expanding the logsinh in eq. (21) as an e-sum, we find a term-by-term agreement with 
eq. (20). It is clear that So = 0, showing that our quantum formulation indeed regularizes the low T anomaly of the 
entropy of ref . [9]. The entropy and specific heat now behave as S ~ C ~ T 3 / 2 , related to the square root singularity 
of the semi-circular law at A = 2J2. This describe gapless excitations. In contrast, the mean field quantum spherical 
ferromagnet does have a gap. [13] 

When also quartet couplings, and possibly higher couplings, are present, replica symmetry is broken. The entropy 
at low T can be obtained by analyzing the continuum limit of the equations for p e . Details of this derivation will not 
be presented here. It turns out that the same behavior occurs, 

Figure 2 presents the specific heat and entropy for a = 1/4 in the situation of Fig. 1, a = 1/2, J2 = 1, J4 = 2. 

In conclusion, we have presented a set of models where the order parameter function can be solved exactly in 
the whole low temperature phase. The critical behavior is exactly the one derived by Parisi for the SK model. As 
compared to that model, the solution is explicit for all T. It is shown that the order parameter function has a universal 
shape in the whole low temperature region. The breakpoint x\ is small near the transition and near T = 0. Well 
below Tq states with very different overlaps, so of very different nature, occur. 

We have also considered a quantum version of the model. It is seen that the shape of the order parameter function 
is not changed. The low T divergency of the entropy in the classical model is eliminated. Now the entropy and specific 
heat vanish as T 3 I 2 . 

It is interesting to point out that the magnetic properties of the model, such as q(x), Xzfc and Xfc have not 
changed qualitatively in the quantum description. The only relevant difference is that the self-overlap qd has become 
a smooth function of T. In our model Xzfc remains non-zero at T = 0. This disagrees with the SK model but is 
very often observed experimentally. Another new aspect of our model is the boundedness of the critical field below 
which freezing takes place. For non-zero field one has in the frozen phase M(H,T) = H / ^/ f"(qo) = M(H) and 
S(T,H) = S(T). Therefore the Parisi- Toulouse hypothesis is satisfied. [19] 

It is hoped that the proposed model can be used as a basis to explain experiments. Hereto inclusion of ferromagnetic 
couplings and extension to vector spins seems of large interest. In a dynamical study of the model a closed set of 
equations will occur that bear the same physics as the more complicated set of the SK-model. 
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FIG. 1. Specific heat and entropy as function of temperature in the classical model. The arrow marks Tcj. 
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FIG. 2. Specific heat and entropy in the quantum model. 
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